jECURITY  CLASSIFICATION  OF  THU  PAGE  fWi  Dma  *nl»r»Q 


|  REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS  I 

BEFORE  COMPLETING  FORM  1 

mnnjumijB 

«■  TITLE  fan*  Submit) 

NEW  RESULTS  ON  PLANAR  TRIANGULATIONS 

>.  TYRE  OF  REPORT  A  PERIOO  COVERED 

Technical  Report  y 

7.  AUTMORfFj 

PETER  D.  GILBERT 

CONTRACT  OR  GRANT  NUMBERfaJ 

NSF  MCS  76-17321  . 
N00014-79-C-0424^ 

S.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Coordinated  Science  Laboratory  / 

University  of /Illinois ^ at  Urbana-Champaign  , 

Ur  bans,  Illinois  61801.  __  > 

10.  PROGRAM  ELEMENT.  PROJECT,  TAM 
AREA  A  WORK  UNIT  NUMBERS 

II.  CONTROLLING  OFFICE  NAME  AND  AOORESS 

/*y Pi  "'  - 

Joint  Services  Electronics  Program 

1 

IS.  NUMBER  OF  PAGES 

42 

14.  MONITORING  AGENCY  NAME  •  AOORESSflf  <11  limn til  trmm  Controlling  Ottlcm) 

IS.  SECURITY  CLASS.  (•!  tfil«  npMJ 

UNCLASSIFIED 

*Sa.  DECL  ASSI  Ft  CATION/ DOWNGRADING 
SCHEDULE 

IS.  OISTRISUTION  STATEMENT  (ol  Ihlm  Htpmt) 

Approved  for  public  release;  distribution  unlimited 

17.  OISTRISUTION  STATEMENT  (o 1  (Ac  mbttrmct  entered  In  Block  20,  It  different  from  Report) 

is.  supplementary  notes 

tl  sj  h) 

-  -  -  .  * 

19.  KEY  WOROS  (Continue  on  reeeree  elde  It  neceeeery  end  Identify  by  block  number) 

Computational  Geometry  Segment  Trees 

Triangulations  Efficient  Algorithms  A/,(  ) 

Greedy  Triangulation  ^  ^  ^ 

20.  ABSTRACT  (Continue  on  rovono  iltfo  If  neceeemry  end  Identity  by  block  number) 

}'  Several  new  results  in  planar  triangulations  are  presented.  In  particular, 

an  GCN^logN)^  time  worst-case  algorithm  for  the  construction  of  the  ’greedy' 
triangulation  is  presented  an  Improvement  over  the  previously  best-known 
construction  in  t ime  >15 (N J).  '  The  minimum  weight  triangulattlon  (MWT)  of  a  simple 
N-vertex  polygon  is  showtT  to  be  constructible  in  time  (KN3) 7/  A  local  property  o' 
MWTs  is  proved,  a  side  result  of  which  is  that  the  shortest  edge  between  a 
pair  of  points  must  be  in  a  MWT.  Also,  it  is  shown  that  the  Delaunay 

00  1473  eoirioNor  i  nov  ss  is  obsolete  UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whm*  Dm •  En(»r*#) 


SIFIED 


SSCgWTV  CLASSIFICATION  OF  TNI*  PAOSflHMa  O «M 


2Q<  ABSTRACT  (continued) 


r Triangulation  is  not  'approximately  optimal',  and  that  the  ratio  of 
'  the  total  length  of  the  Delaunay  Triangulation  of  N  points  to  that  of 
the  MWT  can  be  as  large  as  roughly  N/2. 


Accession  For 

mis 

DOC  TAB 

Urr.".n  ounce  d 
Ju:. tif  ic  !.  1  >•’ 


Dlr.*  rl'  •  v  ■/ 
Av.-ll- 


UNCLASSIFIED 


SSCUWTr  CLASSIFICATION  OF  THIS  NAOtffcn  Omm  bNn« 


NEW  RESULTS  GN  PLANAR  TRIANGULATIONS 


il 


f'q  tr 

sv.  , 


i 


i  (- 


i . 


BY 

.M- 


JPETER  J) ./GILBERT 


U{  '  7  / 


B.S.,  University  of  Illinois,  1976 


THESIS 


Submitted  in  partial  fulfillment  of  the  requirements 
for  the  degree  of  Master  of  Science  in  Computer  Science 
in  the  Graduate  College  of  the 
University  of  Illinois  at  Urb ana -Champaign,  1979 


Thesis  Advisor:  Professor  F.  P.  Preparata 


Urbana,  Illinois 


ABSTRACT 


Several  new  results  In  planar  triangulations  are  presented.  In 
2 

particular,  an  0(N  logN)  time  worst-case  algorithm  for  the  construction 

of  the  'greedy*  triangulation  is  presented,  an  improvement  over  the 

previously  best-known  construction  in  time  O(N^).  The  minimum  weight 

triangulation  (MWT)  of  a  simple  N-vertex  polygon  is  shown  to  be  con- 

3 

structible  in  time  0(N  ).  A  local  property  of  MTTs  is  proved,  a  side 
result  of  which  is  that  the  shortest  edge  between  a  pair  of  points  mist 
be  in  a  MWT.  Also,  it  is  shown  that  the  Delaunay  Triangulation  is  not 
'approximately  optimal',  and  that  the  ratio  of  the  total  length  of  the 
Delaunay  Triangulation  of  N  points  to  that  of  the  MMT  can  be  as  large 
as  roughly  N/2. 
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CHAPTER  1 
INTRODUCTION 

1.1  Computational  Geometry 

Computational  geometry  is  a  relatively  new  field  of  study,  dealing 
with  time  and  space  complexities  of  geometric  problems.  Many  important 
problems,  both  theoretically  and  from  the  point  of  view  of  applications, 
are  geometric  in  nature,  including,  for  example,  wire-routing,  the 
Traveling  Salesman  Problem,  and  the  Steiner  Tree  Problem.  Solving  these 
efficiently  necessitates  developing  new  algorithms  and  new  tools  that 
take  advantage  of  the  geometric  nature  of  the  problems.  It  is  currently 
unclear  to  what  degree  geometry  affects  the  complexity  of  problems,  and 
because  of  the  relative  youth  of  the  subject,  there  remains  much  work 
to  be  done  and  a  large  number  of  open  problems. 

Much  of  the  early  work  in  computational  geometry  has  either  been 
done  or  surveyed  by  M.  I.  Shamos  [1,2],  who  has  given  efficient  algorithms 
for  many  fundamental  problems  in  the  field.  For  a  good  survey  of  work 
in  the  field,  see  reference  [1]. 

This  thesis  is  concerned  with  some  algorithms,  results,  and  techniques 
which  could  prove  very  useful  in  this  rapidly  expanding  area  of  research. 

1.2  Preliminary  Definitions 

Before  proceeding,  it  is  necessary  to  define  some  terms. 

Let  V  be  a  set  of  N  distinct  points  in  the  plane,  and  let  L  be  the 
set  of  edges  (or  line  segments)  between  points  in  V.  Two  distinct 

edges  e^  and  e^  properly  intersect  if  e^  and  e ^  Intersect. at  a  single 
point  which  is  not  an  endpoint  of  both  e^  and  e^. 
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A  triangulation  of  V  is  a  maximal  subset,  T,  of  L  such  that  no  two 
edges  of  T  properly  intersect.  The  term  convex  hull  is  used  to  refer  to 
the  convex  polygon  on  the  boundary  of  a  triangulation  of  a  set  of  points. 
More  precisely,  the  points  on  the  convex  hull  of  a  set  of  N  points  are  all 
those  points  that  are  not  in  the  interior  of  say  triangle  formed  by  any 
three  of  the  N  points.  The  edges  on  the  convex  hull  are  all  the  edges 
with  endpoints  on  the  convex  hull,  and  which  do  not  properly  Intersect  any 
of  the  edges  connecting  pairs  of  points.  Thus,  the  edges  on  the  convex 
hull  must  be  in  any  triangulation  of  the  set  of  points. 

Some  properties  of  triangulations  are  immediately  apparent  or  easily 
derived : 

1.  Each  interior  region  of  a  triangulation  determined  by  T  and 
V  is  triangular  in  shape. 

2.  Each  edge  in  L  is  either  in  T  or  properly  intersects  an 
edge  in  T. 

3.  The  number  of  edges  in  a  triangulation,  |t] ,  is  equal  to 

3N  -  | CH|  -  3,  where  |ch|  is  the  number  of  points  on  the  convex 
hull  of  V. 

Figure  1  should  help  in  the  visualization  of  these  concepts. 


A  triangulation  of  A  convex  hull, 

a  set  of  points 


Figure  1.  N  -  10;  ll|  -20;  |  CH|  -7. 


•'-n 
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1.3  The  Minimum  Weight  Triangulation  Problem 

For  any  given  set  of  points,  there  is  in  general  more  than  one  possible 
way  in  which  they  can  be  triangulated.  One  particular  triangulation  of  a 
set  of  points,  the  Minimum  Height  Triangulation  (also  called  the  Optimum 
Triangulation),  is  defined  as  follows.  Of  all  possible  triangulations  of 
a  set  of  points,  a  triangulation  for  which  the  sum  of  the  lengths  of  the 
edges  is  a  minimum  is  called  a  Minimum  Weight  Triangulation  (MWT)  of  the 
set  of  points. 

Previously,  very  little  has  been  known  about  MfTs  besides  the 
definition,  most  results  being  counter-examples  of  propositions  about 
MWTs,  and  thus  illustrating  what  a  MWT  is  not.  In  Chapter  3  we  present 
a  polynomial- time  algorithm  for  constructing  the  MWT  of  a  simple  polygon, 
and  prove  a  local  property  of  MWTs,  a  side  result  of  which  is  that  the 
shortest  edge  must  be  in  the  MWT. 

E.  L.  Lloyd  [3]  has  conjectured  that  the  problem  of  constructing  the 
MWT  of  a  set  of  N  points  is  HP-complete,  and  has  presented  counter-examples 
to  several  propositions  regarding  MWTs.  Other  work  on  the  MWT  problem  has 
been  done  by  Diippe  and  Gottschalk  [4],  and  Shamos  and  Hoey  [1]. 

1.4  The  Greedy  Triangulation  Problem 

One  algorithm  purported  to  solve  the  MWT  problem  (in  time  polynomial 
on  the  number  of  points)  was  published  by  Duppe  and  Gottschalk  [4]. 
Unfortunately,  the  algorithm  does  not  necessarily  yield  a  MWT  [3],  but  is 
of  interest,  nevertheless.  The  algorithm,  which  defines  the  Greedy  Tri- 
angulation  (GT) ,  is  as  follows: 


1.  Construct  the  (^)  possible  line  segments  between  pairs  of 
points. 

2.  Order  these  line  segments  by  their  lengths. 

3.  Include  the  shortest  of  these  line  segments  in  the  tri¬ 
angulation,  and  eliminate  all  line  segments  which  intersect  this 
shortest  one. 

4.  Continue  constructing  the  triangulation  as  in  step  3,  using 
the  remaining  line  segments. 

An  alternate  description  of  the  construction  of  the  Greedy  Tri¬ 
angulation  is  the  following: 

1.  Set  GT  to  the  empty  set.  Sort  the  edges  of  L  by  length. 

2.  Consider  each  edge  e  in  L  in  order,  from  shortest  to  longest. 

If  e  does  not  properly  intersect  any  edge  already  in  GT,  then  add  e 
to  GT. 

3.  GT  is  the  set  of  edges  in  the  Greedy  Triangulation. 

It  is  clear  that  a  straightforward  implementation  of  this  algorithm 
3  2 

uses  0(N  )  time  and  0(N  )  space.  G.  K.  Manacher  and  A.  L.  Zobrist 

[5]  have  presented  an  O(N^)  time  (worst  case)  and  0(N)-space  algorithm 

2 

to  construct  the  GT,  with  average  running  time  of  0(N  logN).  They  also 

show  that,  contrary  to  what  had  been  thought,  the  GT  need  not  be  a  good 

approximation  to  the  MWT;  that  the  sum  of  the  lengths  of  the  edges  of  the 

1/3 

GT  may  be  greater  than  that  of  the  MWT  by  a  factor  of  0(N  ).  Chapter  2 

2 

presents  two  0(N  )  space  algorithms  for  the  construction  of  the  GT,  with 
worst-case  running  times  of  Offl5'2)  and  0(N^logN),  respectively.  It  is 


also  shown  that  the  total  length  of  the  Delaunay  Triangulation  [2]  can  be 
N/2  times  as  large  as  that  of  the  MWT. 


We  note  also  that  it  has  been  shown  [2]  that  Q(NlogN)  is  a  lower 
bound  on  the  time  required  to  construct  any  triangulation  in  ’good*  form, 
i.e. ,  all  points  which  are  edge-adjacent  to  a  particular  point  are  sorted 
by  polar  angle  around  that  point,  and  that  this  bound  is  achievable. 
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CHAPTER  2 

GREEDY  TRIANGULATIONS 

2.1  Overview 

In  this  chapter  a  "greedy"  algorithm  for  triangulating  a  set  of 

points  is  discussed,  and  two  algorithms  for  the  construction  of  this 

triangulation  are  presented.  Each  of  these  algorithms  is  an  improvement 

3 

over  the  previously  best  known  0(N  )  time  worst-case  construction  of  the 
Greedy  Triangulation. 

A  brute  force  construction  of  the  Greedy  Triangulation  (GT)  can  be 
3  2 

achieved  in  0(N  )  time  and  0(N  )  space  (see  section  1.3).  One  would 
naturally  expect  lower  time  and  space  bounds  for  two  basic  reasons: 

1.  An  0(NlogN)  time  triangulation  algorithm  exists  (for  the 
Delaunay  Triangulation),  this  being  the  lower  bound  for  producing  any 
triangulation  in  'good'  form  (2],  and 

2.  The  GT  problem  is  a  relatively  simple  and  straight-forward 
problem  in  computational  geometry. 

Manacher  and  Zobrist  have  called  triangulations  of  the  following 
form  "conformal"  [5] : 

1.  T  ♦*  t\  L  *-  {all  edges  connecting  pairs  of  points}. 

2.  Determine  an  edge  e€L,  to  be  added  to  the  triangulation t 
and  set  T  «-  T  U  { e} . 

3.  Repeat  step  2  until  T  is  a  triangulation. 

In  a  conformal  triangulation  algorithm,  an  edge  which  is  added  to  the 
triangulation  is  not  deleted  from  the  triangulation  later.  In  this  respect, 
the  two  algorithms  presented  in  this  chapter  are  conformal. 


We  shall  now  present  an  algorithm  for  the  construction  of  a  GT  of 

a  set  of  N  points  in  the  plane.  It  is  similar  in  structure  to  the 

algorithm  of  section  1.3,  but  'latches"  updates  to  the  data  structures 

5/2 

in  order  to  achieve  an  0(N  )  time  worst-case  construction. 

At  the  heart  of  any  algorithm  are  its  data  structures.  In  this 
proposed  algorithm  for  constructing  the  GT  of  N  points,  we  have  L,  which 
is  an  ordered  set  of  the  edges  between  pairs  of  points,  sorted  in 

order  of  increasing  length,  from  which  we  consider  edges  one  at  a  time. 

Also,  we  have  an  NXN  adjacency  matrix  ADJ,  the  elements  of  which  are 
Boolean  values  ADJ(i,j),  initialized  to  1  (iff  ijtj),  and  set  to  0  when 
the  edge  P^Pj  Is  determined  to  not  be  in  the  GT,  so  that  at  the  end  of 
the  algorithm  ADJ  is  a  standard  adjacency  matrix  representation  of  the 
edges  in  the  GT. 

Associated  with  each  point  p  is  a  star-structure,  STAR  ,  which  is 

s  s 

an  ordered  subset  of  V-{p  },  the  points  being  ordered  by  polar  angle 

s 

around  p  .  Initially,  STAR  *  V-{p  },  ordered  by  polar  angle.  If  a  point 
s  s  s 

p^  is  deleted  from  the  star-structures  STARa>  this  corresponds  to  the  deter¬ 
mination  that  the  edge  P8P^  is  not  in  the  GT.  In  the  course  of  the  algorithm, 

the  elements  ADJ(s,i)  which  equal  1  correspond  exactly  with  those  points  p^ 

2 

in  STARg.  Initializing  these  data  structures  can  be  done  in  time  0(N  logN). 

To  achieve  a  "low"  upper  bound  on  the  running  time,  we  "batch"  up¬ 
dates  to  these  data  structures,  so  that  we  process  a(N)  edges  at  a  time. 

These  edges  have  been  determined  to  be  in  the  GT,  and  are  kept  as  an 
unordered  set  of  edges  H,  with  | h|  <  a(N).  In  the  subsequent  analysis,  the 
value  of  a(N)  will  be  chosen  for  optimality. 
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To  Cest  an  edge  p^p^  for  inclusion  in  GT,  we  muse  check  ADJ(i,J), 
which  was  current  as  of  the  last  batched  update,  but  we  oust  also 
check  for  proper  intersection  between  p^p^  and  each  of  the  edges  currently 
in  the  set  H,  which  have  been  added  to  GT  since  the  last  update  of  the 
star-structures  and  the  adjacency  matrix. 

We  are  now  ready  to  sketch  the  algorithm  (see  Figure  2): 

1.  (Initialization) 

Let  V  be  a  set  of  N  points  in  the  plane:  p,,  p_,  p„. 

1  i.  N 

L  *-  {all  edges  whose  endpoints  are  in  v},  sorted  by  length. 

STAR  t-V-Cp  },  sorted  by  polar  angle  around  p  ;  for  s»l,2,...,N. 

3  S  S 

ADJ  (i,  j)  «-  1  for  i,j  -  1,2,...,N,  i*  j.  GT  «-  0;  H  «-  0. 

2.  (Should  a  considered  edge  be  added?) 

Consider  each  edge  p^p^  in  L,  in  order  of  increasing  length. 

If  all  edges  have  already  been  considered,  then  perform  step  4 
(update  data  structures)  a  final  time  and  stop.  If  ADJ(l,j)  -  1, 
and,  if  p^Pj  does  not  properly  intersect  any  edge  in  H,  then 
perform  step  3  (add  p^p j  to  GT  and  H) ,  and  possibly  step  4 
(update  STARs);  otherwise,  continue  considering  edges. 

3.  (Insert  edge) 

GT  ♦-  GT  U  {  PtP  j  3  »  m-HUtp^j}. 

If  | h|  ■  a(N)  then  perform  step  4,  and  continue  with  step  2, 
else  just  continue  with  step  2. 

4.  (Update  Stars  and  ADJ) 

For  all  s€{l,2 . N}  ,  set  STAR  ♦- STAR  -{ p.  |  edge  p  p.  properly 

intersects  an  edge  in  h} . 

If  p  p.  properly  Intersects  an  edge  in  H,  then  set  ADJ(s,i)  to  0 

S  1 

Set  H+*0  and  go  to  step  2. 


A 


d.  A  set  of  3  edges 
In  H. 


e.  STAR,  after  being  f.  STAR...  after  being 

updated  by  H.  updated  by  H. 


Figures  2(a-i).  Examples  of  the  updating  of  STAR- structures 


Most  of  the  work  in  the  algorithm  is  done  in  step  4,  updating  the 
data  structures,  which  we  now  analyze  in  more  detail. 

Given  a  ray  r  from  a  point  p,  and  a  set  E  of  nonintersecting  edges, 
we  let  ”  {e|e€E,  e  and  r  properly  intersect).  Then,  for  two  edges 
el,e2^^r  we  Sa^  el  c^oser  e^  (e^  <  e^)  if  the  intersection  of 

e^  and  r  is  nearer  to  p  than  is  the  intersection  of  e^  and  r.  This  is 
a  proper  ordering,  since,  if  r  and  r1  are  rays  from  point  p  which  both 
intersect  e^  and  then  <  e^  and  e^  <,  e^  would  imply  a  proper 
intersection  between  e^  and  e^ ,  which  we  know  to  be  impossible  by  the 
definition  of  E. 

Lemma  2.2.1:  Let  h  ■«  | h]  .  The  step  "STAR  *-STAR  -  £ p  ]  edge  p  p  properly 

Intersects  an  edge  in  h}"  can  be  done  in  time  0(h  log  h  +  N),  as  can 

the  corresponding  update  of  ADJ. 

Proof;  Consider  the  following  method  of  performing  the  step. 

Sort  the  2h  endpoints  of  edges  in  H  by  polar  angle  around  p  , 

s 

associating  with  each  endpoint  information  saying  which  edge  it  is  from 
(if  any  endpoint  coincides  with  pg,  then  its  corresponding  edge  can  be 
Ignored  in  the  updating  of  this  star-structure,  since  it  will  not  properly 
intersect  any  other  p  p,  ).  Let  H  be  the  sequence  of  edges  from  H  which 

S  K  Tq 

are  intersected  by  a  chosen  ray  r^,  emanating  from  pg,  ordered  according 
to  their  closeness  to  p  ;  H  is  organized  as  a  height-balanced  binary 

S  Tq 

tree  (AVL  priority  queue).  Initially  we  set  U«-H  ,  where  rn  is  a  ray 

r0  U 

from  pg  at  polar  angle  0.  Then  we  scan  the  points  of  HUSTARg  sorted 
according  to  their  polar  angle  around  point  pg. 

Let  p^  be  the  point  being  scanned.  Either  p^€  STARg  (Fig.  3a, b)  or 
Pk  is  an  endpoint  of  one  or  more  edges  in  H,  (Fig.  3c,d,e)  or  both.  If 
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a.  Edge  p  pfc  properly  inter¬ 
sects  the  closest  edge  in 
U  (solid  edges  are  in  U). 


b.  Edge  p  p,  does  not  properly 

intersect  the  closest  edge  in  U 
(solid  edges  are  in  U). 


c.  Point  p,  is  a  'leading' 
endpoint  of  an  edge  in  H. 


d.  Point  p,  is  a  'trailing' 
endpoint  of  an  edge  in  H. 


e.  Point  p,  is  both  a  'leading' 
and  a  'trailing'  endpoint  of 
two  edges  in  H. 


Figure  3(a-e).  Scanning  points  from  STARS  and  endpoints  of  edges  in  H 
(see  lenraa  2.2.1). 


Pk£STARg  and  pgpk  properly  intersects  the  closest  edge  in  U  (Fig.  3a), 


then  set  STARs  «-  STARg  -  {p^}  (and  ADJ(s.k)  is  set  to  0). 

Additionally,  if  pk  is  an  endpoint  of  an  edge  in  H,  then  U  oust  be 

updated.  If  pk  is  the  'leading*  endpoint  (Fig.  3c),  in  the  given  ordering 

of  polar  angles,  of  an  edge  p.p.  in  H,  then  p  p .  must  be  added  to  U;  if  p, 

K  J  K  J  lC 

is  the  'trailing'  endpoint  (Fig.  3d),  then  PjPk  must  be  removed  from  U. 

Since  p,  is  a  trailing  endpoint  of  the  edge  p.p.  in  H  if  and  only  if  p.p. 

K  J  K  J  K 

is  already  in  U,  this  update  to  U  can  be  written: 

U  «-  UU{PjPk|PjPk€  H,  PjPk?U} 

pjPk6u) 

We  continue  scanning  until  all  points  in  STARg  have  been  scanned.  Since  an 
edge  is  added  to  and  deleted  from  U  only  once,  the  time  bound  follows 
immediately.  □ 

Notice  that  with  this  GT  algorithm,  at  the  algorithm's  completion,  the 
Greedy  Triangulation  is  available  in  three  fortrs:  (1)  As  a  set  of  edges  (GT) : 
(2)  as  an  adjacency  matrix  (ADJ) ;  (3)  as  a  collection  of  star-structures,  where 

each  point's  star-structure  contains  that  point's  neighbors,  in  polar  order. 
Analysis  of  GT  algorithm: 

Initializing  L  (sorting  ^N(N-l)  lengths):  0(N  logN)  time. 

Initializing  the  N  star-structures:  0(N*NlogN)  time. 

2 

Initializing  the  adjacency  matrix:  0(N  )  time. 

Checking  an  edge  for  inclusion  in  GT  (checking  the  edge  against 
all  current  members  of  H) :  0(a(N))  time. 

Checking  ^N(N-l)  edges:  0(N^a(N))  time. 


v* 
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Updating  a  star-structure  with  a(N)  edges  (by  lemma  2.2.1): 

0(a (N) log  a(N)  +  N)  time. 

Performing  0(3N/a(N))  batched  updates  (a  triangulation  contains 

no  more  than  3N  edges),  of  N  star-structures  each: 

2 

0(N  (a(N)log  a(N)  +  N)/a(N))  time. 

Overall:  0(N2logN  +  N2a(N)  +  N3/a(N))  time. 

By  choosing  a(N)  *  O(N^),  where  the  constant  of  proportionality 
is  independent  of  N,  we  obtain  the  following: 

Theorem  2.2.2 :  The  Greedy  Triangulation  of  a  set  of  N  points 
can  be  constructed  in  (KN572)  time  and  0(N2)  space. 

2 

2.3  An  Q(N  loglO-time  Greedy  Triangulation  Algorithm 

In  this  section  we  describe  another  algorithm,  similar  in  structure 

to  that  of  section  2.2,  for  constructing  the  GT  of  a  set  of  N  points. 

A  new  data  structure  is  used  which  will  permit  checking  an  edge  for 

inclusion  in  the  GT  in  time  O(logN)  and  the  updating  for  a  single  edge 

2 

added  to  the  GT  in  time  O(NlogN),  thereby  achieving  an  0(N  logN)  running 
time  for  the  algorithm. 

The  data  structure  is  an  array  of  Range-trees.  For  Integers 
a,b€ { 1,2, . . . ,n}  a  Range  tree  of  the  range  (a,b)  (RT(a,b))  is  a  binary 
tree  whose  nodes  contain  a  subrange  (of  (a,b))  field,  such  that: 

1.  the  root  has  subrange  (a,b). 

2.  if  the  subrange  in  a  node  has  width  1  (a+l»b),  then  the  node 
has  no  subtrees. 

3.  otherwise,  the  node  has  left  and  right  subtrees,  which  are 
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-  ■  .  ier-  -,-  .  .. 


RT(a,  L(a+b)/2J)  and  RT (L(a  +  b)/2j  ,  b),  respectively. 

It  is  simple  to  show  that  RT(a,b):  is  a  binary  tree,  each  node 
having  0  or  2  sons;  has  (b-a)  leaves;  has  2(b-a)-l  nodes  (inner  &  leaves); 
is  balanced;  and  has  maximum  depth  of  Tlog^Cb-a)!  +1  (see  Figure  4).  We 
call  basis  ranges  the  subranges  specified  in  the  nodes  of  a  Range-tree. 

An  exact  covering  of  a  range  (aQ.bg)  is  a  set  of  basis  ranges  (a^.b^)  such 
that  the  union  of  the  basis  ranges  is  equal  to  (ag.bg). 


Figure  4.  A  Range -tree,  RT(1,12),  with  the  minimal  exact  cover 
of  the  range  4.. 11  marked. 


Claim:  Any  subrange  (ag.bg)  of  the  range  (a.b)  can  be  exactly  covered 
using  no  more  than  2Tlog2(b-a)l  +  2  basis  ranges  from  RT(a.b). 

Proof:  Consider  all  basis  ranges  (a^b^)  of  RI(a,b)  such  that 
aQ  <  a^  <  b^  <  bg.  Let  this  be  the  initial  covering,  which  we  now  prune 


E 

[ 

S 

l 


to  the  claimed  number  of  basis  ranges 


All  of  Che  basis  range  nodes  which  are  leaves  of  Che  Range-cree 
have  brothers  which  are  also  in  Che  initial  covering,  except  possibly 
two  such  nodes  at  the  extremes  of  the  range.  A  pair  of  brothers  is 
exactly  covered  by  their  father.  Remove  these  pairs  of  brothers  from 
the  covering  set,  leaving  at  most  two  basis  ranges  at  this  level. 
Proceeding  similarly  level-by-level  to  the  root,  it  is  clear  that  any 
subrange  of  (a,b)  can  be  exactly  covered  using  basis  ranges  in  number  no 
more  than  twice  the  depth  of  the  Range-tree,  which  completes  the  proof.  [ 
Also,  it  should  be  realized  that  all  these  nodes  can  be  marked  or 
visited  in  time  O(log(b-a)). 

Each  ps  will  have  associated  with  it  a  RT(1,N),  the  nodes  of  which 

also  contain  a  line-segment  field,  described  below.  A  number  k  as  the 

end  of  a  basis  range  in  this  Range-tree  corresponds  to  the  direction  from 

p  to  another  point.  If  the  N-l  points  V-{p  }  are  sorted  by  polar  angle 
s  s 

around  p  ,  and  p,  is  the  k1"*1  point  in  this  ordering,  then  the  number  k  as 
S  1 

the  end  of  a  basis  range  is  associated  with  point  p^  in  the  numbering 

scheme  around  p  .  An  N X  N  array,  Global-to- local,  is  used  to  convert 
s 

indices  from  the  global  numbering  to  the  local  numbering  scheme  around  pg 

In  the  above  case,  we  have  Global-to-local  (s,i)*k.  -This  array  can  be 
2 

initialized  in  0(N  logN)  time  (N  sorting  operations,  of  N-l  points  each). 

Direction  N  (around  p  )  is  considered  to  be  the  same  as  direction  1,  for 

s 

circularity. 

The  line-segment  field  of  a  node  with  basis  range  (a^,b^),  if  not 
empty,  will  contain  the  name  of  the  edge  which  has  been  added  to  the  GT, 
spans  the  range,  and  is  the  closest  to  p  .  Thus,  suppose  that  edge  p  p 
is  being  considered  for  inclusion  in  the  (7T;  we  first  refer  to  the 
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Global-to-local  array  and  obtain  an  integer  k  ■  Global-to- local  (s,i), 

and  then  traverse  a  path  of  the  Range-tree  for  point  p  from  the  root 

to  a  leaf  corresponding  to  the  integer  k.  If  p  p  properly  Intersects  any 

8  X 

of  the  edges  in  the  line-segment  field  of  the  nodes  of  this  path,  then 
p^p^  cannot  be  added  to  the  GT:  otherwise  (since  edges  are  considered  in 
order  of  increasing  length),  the  edge  p  p,  should  be  added  to  the  GT. 

Because  the  depth  of  the  tree  is  O(logN),  testing  an  edge  for  inclusion 
in  GT  requires  O(logN)  time. 

When  an  edge  is  determined  to  be  in  the  GT,  the  line-segment  fields 

of  the  Range-trees  must  be  updated.  Consider  updating  the  Range-tree  for 

point  p  with  the  edge  p  p  .  If  i  ”  s  or  j*s,  then  no  update  of  this  RT 
s  i  j 

need  be  done.  Otherwise,  we  look  up  Global-to-local  (s,i)  and  Global-to- 

local  (s,j),  and  let  f  and  g  be  the  smaller  and  larger  of  these,  respectively. 

Then,  visiting  each  node  of  the  RT  in  the  minimal  exact  cover  of  (f,g),  if 

the  line-segment  field  is  empty,  or  p  p  is  closer  to  p  than  the  edge 

l  J  s 

named  in  that  field,  set  the  line-segment  field  to  P^j  (actually,  if  a 

ray  from  p  at  polar  angle  0  properly  intersects  p.p.,  then  we  must  break 
s  1  J 

(f,g)  into  the  two  parts  (l,f)  and  (g,N),  and  perform  a  similar  update  with 

each  of  these).  This  can  be  done  in  time  O(logN). 

To  analyze  this  proposed  GT  algorithm,  we  note  that  the  initialization 

of  L,  the  N  Range-trees,  and  the  array  Global-to-local  can  be  done  in 
2 

0(N  logN)  time.  Checking  an  edge  for  inclusion  in  GT  takes  O(logN)  time, 
and  hence,  checking  ^N(N-l)  edges  for  inclusion  in  GT  takes  0(N  logN)  time. 
Since  updating  a  single  Range-tree  with  a  single  edge  can  be  done  in  time 
O(logN),  updating  N  Range-trees  with  the  0(N)  edges  added  to  the  GT  can  be 
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done  in  O(NlogN)  time.  Thus,  the  overall  running  time  of  the  algorithm 
2 

is  0(N  logN),  and  we  have  proved: 

Theorem  2.3:  The  GT  of  a  set  of  N  points  can  be  constructed  in 

2  2 
0(N  logN)  time  and  0(N  )  space. 

By  being  able  to  essentially  occlude  (possibly)  many  points  from  being 

visible  from  a  point  p  at  once,  without  sacrificing  fast  checking  of  edges 

s 

for  inclusion  in  the  GT,  we  achieve  a  "nicely  balanced"  O(N^logN)  time 
construction  of  the  Greedy  Triangulation  for  N  points. 

We  note  here  that  a  Range-tree  structure  has  been  useful  in  several 
other  computational  geometry  algorithms  [8,9],  such  as  computing  the  area 
of  the  union  of  N  aligned  rectangles  (sides  parallel  to  the  x  and  y  axes) 
in  time  O(NlogN)  [10]. 
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CHAPTER  3 

MINIMUM  WEIGHT  TRIANGULATIONS 

3. 1  Introduction 

Presented  in  this  chapter  are  several  results  concerning  Minimum 
Weight  Trianguiaticns  (MWT). 

A  construction  of  the  MWT  of  a  simple  polygon  in  O(N^)  time  and  O(N^) 
space  is  demonstrated  in  section  3.2. 

The  major  result  of  this  chapter  involves  a  proof  of  i  local  property 
of  MWTs,  which  may  be  very  useful  in  further  work  on  MWTs.  A  side  result 
of  this  theorem  answers  the  previously  open  question  of  whether  the  shortest 
edge  is  in  a  MWT  of  a  set  of  vertices;  it  is.  We  now  informally  describe 
this  theorem. 

Given  two  simple  closed  curves  in  the  plane,  one  entirely  enclosing 
the  other,  the  plane  is  divided  into  three  regions;  inner,  middle,  and  outer. 
We  consider  MWTs  of  a  set  of  points  which  lie  in  the  inner  and  outer  regions 
only;  that  is,  there  are  no  points  in  the  middle  region.  Intuitively,  if 


Figure  6.  Illustration  of  a  separator  segment. 


21 


Inner  region  is  correctly  placed,  and  is  sufficiently  "small",  we  would 
expect  there  to  be  no  edge  in  the  MWT  connecting  two  points  in  the  outer 
region,  passing  through  the  inner  region,  and  separating  the  points  in 
the  inner  region  into  two  nonempty  sets.  We  call  such  an  edge  a  separator 
segment. 

The  theorem  is  primarily  concerned  with  the  case  when  the  two  curves 
are  a  pair  of  concentric  circles,  and  on  sufficient  conditions  which  guar¬ 
antee  that  there  be  no  separator  segment  in  the  MWT  of  these  points.  The 
shortest  edge  between  a  pair  of  points  is  easily  shown  to  meet  these  require¬ 
ments,  from  which  the  fact  that  this  shortest  edge  must  be  in  a  MWT  of  the 
points  follows. 

3.2  The  MWT  of  a  Polygon 

We  state  the  problem  of  triangulating  a  polygon  as  follows:  "Given 
an  N-vertex  simple  polygon  P,  subdivide  its  interior  into  triangles  whose 
vertices  are  also  vertices  of  the  polygon."  Correspondingly,  of  all  possible 
triangulations  of  a  simple  polygon,  a  triangulation  for  which  the  sum  of  the 
lengths  of  the  edges  is  a  minimum  is  called  a  Minimum  Weigvu  Triangulation 
of  the  polygon.  This  problem  is  simpler  than  that  of  constructing  a  MWT  of 
an  arbitrary  set  of  points  because  there  are  no  points  in  the  interior  of 
the  region  to  be  triangulated. 

For  convenience,  we  consider  at  first  the  problem  of  constructing  the 
MWT  of  the  vertices  of  a  convex  polygon.  For  simplicity,  we  call  such  sets 
of  points  void  hulls  (Figure  7). 
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Figure  7.  A  void  hull  of  10  points. 

Proposition:  The  MWT  of  a  void  hull  of  N  points  can  be  constructed 

3  2 

in  0(N  )-time  and  0(N  )-space. 

Proof :  The  proof  is  algorithmic  and  is  a  dynamic  programming  argument. 
We  label  the  points  Pq,p^,...,pn  ^ ,  ordered  clockwise  around  the  void  hull, 
and  consider  the  MWT  of  the  s+1  points  p^ »Pi+1 > • • • »P^+S »  where  the  addition 
in  the  subscripts  is  modulo  N.  The  edge  p^p^+g  must  be  in  the  MWT  of  these 
points,  since  it  is  on  the  convex  hull  of  these  points,  and  if  s+1  is  no  less 
than  3,  then  P^P^+S  must  be  a  side  of  a  triangle,  the  third  vertex  of  which 
must  be  one  of  the  points  P^+^ »Pi+2 ’ *  *  * ’ Pi+s - 1  (Fi8ure  8). 


Figure  8. 


Pi+j 


Illustration  of  the  dynamic  programming  argument 
for  the  MWT  of  a  void  hull. 


Assume  that  the  third  vertex  of  this  triangle  in  the  MWT  is  P1+j>  where 
j  may  range  from  1  to  s-1.  By  the  optimality  principle  of  dynamic  programing 

[7],  the  triangulations  of  the  j+1  points  P£>Pi+2 . ^i+j  ’  t*ie  8+^“j 

points  P^+j > *  *  * »P^+S  must  Also  be  of  minimum  weight.  Thus,  if  the  third 
vertex  of  this  triangle  is  Pi+j»  then  the  length  of  the  MWT  of  the  s+1  point 
is  the  sum  of 

the  length  of  the  MWT  of  the  j+1  points  P^»P^+i » • • • >P^+j » 
the  length  of  the  MWT  of  the  s+l-j  points  P^+j > • • • »P^+S 
and  the  length  of  the  edge  P^>Pi+8* 

We  notice  that  each  of  these  MWT  problems  is  smaller  than  the  (s+l)-point 

problem.  By  considering  each  of  j -  1,2 . s-i  in  turn,  calculating  the  MWT 

length  under  each  such  value  of  j,  and  taking  the  minimum  of  these,  the  MWT 

problem  for  the  s+1  points  P^'P^^ » •  •  •  »P^+S  is  reduced  to  smaller  MWT  problems 

Reconstruction  of  the  MWT  is  facilitated  by  storing  along  with  the  length  of 

this  MWT,  the  value  of  j  for  which  the  minimum  is  achieved. 

By  solving  all  the  2-consecutive  point  MWT  subproblems  (s-1),  then  all 

the  3-consecutive  point  MWT  subproblems  (s-2) . and  finally  the  N-con» 

3  2 

secutive  point  MWT  problem,  we  obtain  an  0(N  )-time  and  0(N  ) -space  algorithm 
for  the  construction  of  the  MWT  of  a  void  hull  of  N  points. 

It  should  be  clear  that  this  method  is  easily  extended  to  any  simple 
polygonal  region;  when  constructing  the  MWT  of  the  s+1  points  p^,pi+^, . . . »Pi+g 
if  the  edge  P^P^+8  i-3  not  completely  within  the  polygon,  then  we  simply  set 
the  MWT  length  of  this  (s+l)-polnt  subproblem  to  Infinity.  This  augmented 
algorithm  has  the  same  time  and  space  bounds. 

Unfortunately  no  dynamic  programing  approach  has  so  far  been  found  to 


yield  a  polynomial  algorithm  for  more  general  MWTs. 


We  note  in  passing  that  this  discussion  makes  clear  the  following 
recurrence  on  the  number  of  ways  to  triangulate  a  void  hull  of  N  points 

(P^);  for  any  k  (2  5^11-1),  a  void  hull  of  N-points  can  be  triangulated 

in  a  number  of  ways  which  is  the  product  of  the  numbers  of  trianguiations 
of  void  hulls  with  k  and  N+l-k  points,  respectively;  thus: 

PN  ‘  P2PN-1  +  P3PN-2  +  •••  +  Pn-1P2’  11  >  2|  and 

P2  •  '■ 

The  solution  of  this  recurrence  [6]  is 


2N-4 


N, 


N  “  (  N-2  }  1  (N‘l)  *  4 /16N(ttN) 


1/2 


The  values  of  P^  through  P^  are  1,2,5,14,42,  the  well  known  sequence  of 
Catalan  numbers. 


3.3  Annular  Decision  Criteria 

Consider  a  set  of  points  in  the  plane,  with  the  property  that  a  subset 
of  these  points  is  enclosed  in  two  concentric  circles  with  radii  r^  and  r^ , 
where  r^  <  r^ ,  and  there  are  no  points  in  the  annulus  formed  by  the  two 
circles.  The  annulus  divides  the  plane  into  three  regions  (inner,  middle 
and  outer),  and  no  points  are  in  the  middle  region.  Suppose  that  in  the 
MWT  of  the  set  of  points  there  is  an  edge  connecting  a  pair  of  points  in  the 
outer  region,  passing  through  the  inner  region,  and  dividing  the  points  in 
the  inner  region  into  two  non-empty  sets.  As  mentioned  in  the  introduction 
such  an  edge  is  called  a  separator  segment  (SS).  If  the  outer  radius  (r^) 
were  sufficiently  large  as  compared  with  the  inner  radius  (r^),  one  might 
intuitively  expect  that  there  be  no  separator  segment  in  the  MWT  of  the 
full  set  of  points  (Figure  9). 


Figure  9.  Example  of  a  separator  segment. 

We  now  consider  the  question  of  how  large  the  ratio  of  the  radii  must 
be  in  order  to  guarantee  that  there  by  no  separator  segment  in  the  MWT. 

We  assume  without  loss  of  generality  that  ■  1  and  let  r^/r^  “  r' 

Our  approach  is  by  contradition.  We  assume  that  a  given  triangulation 
be  a  MWT  with  separator  segments,  and  then  show  by  construction  that  there 
is  a  shorter  length  triangulation  with  no  separator  segments.  Specifically, 
the  separator  segments  are  removed,  thereby  obtaining  some  untriangulated 
polygonal  regions;  next,  these  polygonal  regions  are  re-triangulated  by 
inserting  edges  whose  total  length  is  smaller  than  that  of  the  removed 
separator  segments  (Figure  10). 


Example  of  Separator  segments  A  shorter  length  tri¬ 
separator  segments  removed  angulation  with  no 

separator  segments 

Figure  10.  Example  of  Re-triangulation. 

It  should  be  clear  that  each  of  these  polygonal  regions  may  be  re¬ 
triangulated  individually,  and  that  each  of  them  contains  exactly  two 
vertices  in  the  inner  region,  which  we  call  p  and  p  .  Calling  p  and 

A  o  A*  l 

PA+^  the  vertices  adjacent  to  p^  on  the  polygon,  we  also  note  that  the  edge 

PA-1PA+1  WflS  111  the  alle8ed  MWT  (i-e.  PA_jPA+^  was  an  SS). 

Some  preliminary  lemmas  are  required. 

Lemma  3.4,1:  If  r  ^  Jl  then  the  length  of  any  separator  segment  is 
greater  than  the  length  of  p  p  . 

Proof:  Clearly,  length  p  p  <  2,  and  length  (SS)  >  2,/r^-l  >  2  (Figure  11). 

A  D 


Lemma  3,4.2:  If  r  >  5/3,  then  there  exists  a  polygon  vertex  p  outside  the 

v 

annulus  such  that  for  any  separator  segment,  we  have  length  (SS)  >  length 

(PAP„)»  length  (SS)  >  length  (p_p_) ;  moreover,  each  of  the  segments  p  p 

and  p  p  is  either  an  edge  of  the  polygon  or  lies  in  its  interior. 

D  C 

Proof:  We  select  pc  as  a  polygon  vertex  outside  the  annulus  for  which  the 
distance  to  the  center  of  the  annulus  is  a  minimum.  Let  this  distance  be 


d  (Figure  12). 

From  Figure  12  we  see  that  length  (SS)  >  for  any  separator 

segment,  and  that  by  the  triangle  inequality,  length  (PaPq)  <  d+1  and 

length  (p  p  )  <  d  +  1.  By  simple  algebra,  d  >  r  >5/3  implies  2v/l^-l  >  d  +  1, 
B  C 

from  which  we  have  the  two  inequalities. 

The  argument  that  p  p.  is  either  an  edge  of  the  polygon  or  lies  in  its 

interior  is  also  by  contradiction.  Assume  that  segment  PAP^  either  lies 

outside  the  polygon,  or  properly  intersects  it,  and  consider  the  path  p 

from  p  to  p„  formed  by  edges  of  the  polygon,  and  not  passing  through  p  . 
AC  d 
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Figure  12.  Selecting  p  . 

w 

Since  the  vertices  p  and  p  lie  within  the  circle  of  radius  d,  and  p 

A  o  V-* 

lies  on  this  circle,  the  triangle  p.P-P,,  is  contained  in  this  circle.  Re- 
call  that  any  vertex  of  PAC>  other  than  pA>  and  pc>  must  lie  outside  the 
circle  of  radius  r.  If  any  vertex  on  path  PA(,  lies  within  this  triangle,  it 
must  also  lie  within  the  circle  of  radius  d,  and  its  distance  to  the  center 
of  the  annulus  is  less  than  d.  This  is  a  contradiction,  because  p^  was 
chosen  as  the  vertex  outside  the  annulus  and  closest  to  the  center  of  the 
annulus,  and  it  is  at  distance  d.  Thus  we  are  left  with  two  cases  which 
are  typified  by  the  following  examples. 
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Segment  p^p^  lies  outside  Segment  P^Pq  Intersects 

the  polygon.  polygon. 


Figure  13.  Paths  from  pA  to  pc< 


If  p  p  lies  outside  the  polygon  then  the  first  edge  of  the  path  from 
A  L 

p  to  p  on  the  polygon  must  cross  segment  p  p  .  If  P  P„  intersects  the 
AC  o  L  A  L 

polygon,  then  the  edge  it  intersects  also  crosses  segment  PgPc»  since  no 

vertices  are  within  triangle  PAP_P«.  In  either  case,  since  p  is  an  end- 

ABC  C 

point  of  some  separator  segment,  this  separator  segment  must  properly 

intersect  the  polygon,  but  this  is  also  a  contradiction. 

Therefore,  segment  p  p  must  be  an  edge  of  the  polygon,  or  lie  within 

A  C 

the  polygon,  and  similarly  for  segment  p  p  .  3 

o  L 

Lemma  3.4.3:  Given  a  polygon  P  -  (P^  1P2 » •  •  •  >Pn) »  ant*  a  set  T  of  the  n-3 
interior  edges  of  a  triangulation  of  P,  and  three  distinguished  vertices 
of  the  polygon,  pfl,  p^,  and  p£ ,  such  that  neither  pa  nor  p^  is  an  endpoint 
of  any  edge  in  T,  and  p£  is  any  other  vertex,  there  is  a  one-to-one  onto 
mapping  of  the  edges  in  T  to  the  points  P- { pfl , p^ , pc }  such  that  each  edge 
mapped  to  one  of  its  endpoints  (Figure  14). 
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Pr 

Figure  14.  Mapping  of  edges  to  vertices 

(mapping  is  denoted  by  arrows). 

Proof :  (By  induction).  The  lemma  is  trivial  for  n*4  (triangulation  of 

a  quadralateral) .  Thus,  assume  the  lemma  valid  for  all  k  <  n.  Let  t^  be 
the  number  of  edges  in  T  which  have  as  an  endpoint.  Then 

i 

Z  t  -  2|t|  -  2n-6,  t  - 1  -  0, 

p.€p  1  a  b 

i 

and,  by  the  pigeon-hole  property,  at  least  two  of  the  n-2  points 

P-{p  »Pl5  must  have  t  »  1.  Choose  one  of  these  which  is  not  p  ,  call  it 

H  D  x  C 

p^,  and  let  PrPg  be  the  unique  edge  of  T  incident  to  pf.  Edge  PrPg  Is 
mapped  to  p^. 

Consider  now  the  polygon  p',  whose  vertex  set  is  P-{pr;:  obviously 
T'  =  T-{prPg}  is  a  triangulation  of  P1.  Since  jp'i  »  n-1  <  n,  then  by 
the  inductive  hypothesis,  with  the  same  distinguished  vertices  p^,  p^,  and 
Pc  as  before,  there  is  a  one-to-one  onto  mapping  of  the  edges  of  T'  to  the 
vertices  P * -{ pa , p^ , pc } ,  with  each  edge  being  mapped  to  one  of  its  endpoints, 
and  the  lemma  is  proved,  -j 

We  can  now  describe  the  re-triangulation  of  a  polygon  P.  Assuming 
r  >  5/3,  we  find  a  vertex  p  as  in  lemma  3.4.2,  and  construct  a  one-to-one 

w 

correspondence  of  the  interior  edges  of  the  alleged  MWT  to  the  vertices 
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P-{p, »P„»Po3 ,  as  in  lemma  3.4.3.  Thus,  any  separator  segment  (interior 
A  o  C 

edge)  is  longer  than  each  of  PAPg>  PgPc>  and  PAPC>  lennia  3.4.2,  and 

these  three  edges  divide  the  original  polygonal  region  into  three  smaller 

polygonal  regions  and  the  triangle  P.P„P„. 

A  D  L 


Figure  13.  Allocation  of  edges  and  division  of  polygon. 


Consider,  for  example,  re-triangulating  the  smaller  polygonal  region 
which  has  PAPC  as  an  edge  (the  other  two  polygonal  regions  may  be  handled 
similarly),  and  let  t  be  the  number  of  vertices  of  this  polygon.  If  this 


polygon  is  degenerate,  i.e.  it  consists  of  the  single  edge  p  p  ,  then  no 

A  v 

re-triangulation  occurs:  assume  therefore  that  the  polygon  be  nondegenerute 

and  consider  three  consecutive  vertices  p.  ,,  p.,  p.  .  on  this  smaller 

rJ-l  j  j+1 
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polygon  such  tha  p.  £  ^Pa'P^*  Then  p.  must  have  an  assigned  edge.  If 

the  length  of  the  segment  p^  ^Pj+i  *-ess  than  the  length  of  the  edge 

assigned  to  p  ^ ,  and  p  j.p  ^  lies  within  this  smaller  polygon,  then 

Pj  ^Pj+^  may  he  added  to  the  new  triangulation,  vertex  p^  and  its  associated 

edge  may  be  removed  from  further  consideration,  and  the  process  may  be 

continued.  As  the  final  step,  edge  p  d„,  which  is  shorter  than  any 

A  C 

separator  segment,  is  added  to  the  new  triangulation.  The  problem  is  to 
find  such  points  p^. 


Figure  16.  The  separator  segment  assigned  to  p  is  replaced 
by  a  shorter  edge. 

We  now  show  that  if  r  >  2,  then  a  vertex  of  the  original  polygon  tor 
which  the  distance  d  to  the  center  of  the  annulus  is  a  maximum  may  be 
chosen  as  a  suitable  p^. 

The  edge  p^  ^Pj+i  lies  within  the  polygon,  since  the  vertices  Pj_^  and 
p,  .  must  lie  within  the  circle  of  radius  d  and  centered  at  the  center  of 

J  +  i 

the  annulus,  and  the  angle  P j _ i? j*3 j+i.  ^ess  t'ian  180°  •  We  now  show  that 
the  length  of  p^  ^Pj  +  ^  ^ess  than  the  length  ot  the  edge  assigned  to  p^. 


vn#'" 
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The  edge  assigned  to  print  must  pass  through  the  innermost  circle 
(of  radius  1),  thus  its  length  must  be  greater  than  or  equal  to  J?Ti  +  73 
(Figure  17).  By  our  choice  of  p  ^  ,  points  and  must  lie  within 

the  circle  of  radius  d  and  centered  at  the  center  of  the  annulus.  We  now 
distinguish  three  cases. 

(1)  If  both  these  points  lies  within  the  circle  of  radius  1,  then 


length  (Pj_jPj+1)  <2<2j2<  7d  -1  +  73  ,  since  d  >  2. 

(2)  If  one  of  these  points  lies  within  the  circle  of  radius  1,  then 

by  the  triangle  inequality  length  (p^  ^Pj+]7  —  d  +  l>  an<*  ^  2  implies  that 

d  +  1  <  7d2-l  +  73  . 

(3)  If  both  Pj  ^  and  Pj+^  lie  outside  the  circle  of  radius  1,  they 

must  also  lie  outside  the  circle  of  radius  2  (because  the  annulus  is  by 

hypothesis  empty).  In  Figure  17  we  construct  the  line  Pj_^q  tangent  to 

the  circle  of  radius  1,  with  p.  and  the  center  of  the  annulus  on  the  same 

J 

side  of  this  line,  and  with  q  on  the  circle  of  radius  2  so  that  the  tangency 
point  t  is  between  Pj_^  and  q.  The  separator  segment  assigned  to  the  other 
endpoint  of  ^  must  lie  within  the  half-plane  delimited  by  this  line  which 
includes  the  center  of  the  annulus.  Constructing  Che  other  tangent  through 
q  to  the  circle  of  radius  1,  we  claim  that  Pj+1  must  lie  within  the  wedge 
formed  by  these  two  tangents  and  including  the  center  of  the  annulus 
(Fig.  18),  indeed,  by  constructing  the  tangent  Pj+]/l'>  ®s  p^_^q  was  c°n_ 
structed,  the  separator  segments  assigned  to  Pj.^  and  Pj+i  wou^  properly 
intersect,  which  is  impossible.  Thus  the  largest  value  of  the  distance 

between  p  ^  and  p  ^  occurs  when  they  are  on  the  circle  with  radius  d, 

i.e.  when  Pj.^*ci  and  Pj  +  1  *  C2*  noting  that  the  angle  p^_^q  Pj+1  is 


edges  assigned  to  Pj_^  and  Pj+i  properly  inter 
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60°  and  the  triangle  tj^  q  t^  is  equilateral,  length  <  Vd  -1  + 

■J3  <  length  (separator  segment  assigned  to  Pj). 

Thus,  regardless  of  the  positions  of  vertices  Pj_^  and  pj+^»  length 
(pj_lPj+l> < length  (separator  segment  assigned  to  p^),  and  the  vertex  on 
the  polygon  farthest  from  the  center  of  the  annulus  is  a  suitable  pj. 

We  include  edge  Pj_^Pj+^  In  the  new  triangulation,  no  longer  consider 
vertex  pj  and  its  associated  separator  segment,  and  find  another  suitable 
point  (the  next-farthest  vertex).  Continuing  this  process,  the  original 
polygon  is  re-triangulated  using  a  shorter  total  length,  which  implies 
that  the  alleged  triangulation  was  not  a  MWT,  and  we  may  state  the  following 
theorems. 


Given  a  set  of  points  in  the  plane  V,  and  a  subset  S  of  V  with  the 
properties  that  the  points  of  S  may  be  enclosed  by  a  circle  Cj^  of 
radius  r^  and  the  points  of  V-S  are  outside  the  circle  concentric  with 
and  with  radius  r^  >  r^,  we  have: 

Theorem  3.4.4:  If  r^/r^  2,  then  in  the  MWT  of  V  there  is  no  separator 

segment  connecting  vertices  of  V-S,  passing  through  the  inner  circle,  and 
dividing  the  vertices  of  S  into  two  nonempty  sets. 

If  two  of  the  vertices  in  the  inner  region  are  diametrically  opposite 


on  the  inner  circle,  then  a  similar  analysis  can  be  done,  with  the  ratio 
of  radii  reduced  to  «/3.  Figure  19  corresponds  to  Figure  17  of  the 
previous  analysis.  Figure  20  illustrates  the  case  in  which  the  separator 
segment  does  not  pass  through  these  two  vertices. 


Figure  20.  Separator  segments  not  passing  between 

"A  ,nd  V 


i 
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Theorem  3.4.5:  If  ^/r^  —  an<i  two  vertices  °f  S  are  on  a  diameter  of 
the  inner  circle,  and  on  the  inner  circle  then  in  the  MWT  of  V  there  is 
no  separator  segment  connecting  vertices  of  V-S,  passing  through  the  inner 
circle  and  dividing  the  vertices  of  S  into  two  nonempty  sets.  XXX 

As  a  side-result,  we  have: 

Corollary  3.4.6:  The  shortest  edge  is  in  a  MWT  of  any  set  of  points. 

Proof:  Let  pgpt  be  the  shortest  edge,  and  let  2d  be  its  length,  as  in 

Figure  21.  Since  pgpt  is  shortest,  this  implies  that  there  are  no  other 

vertices  in  the  circle  centered  at  pg  with  radius  2d,  and  similarly  for 

the  circle  of  radius  2d  centered  at  p^.  The  annulus  of  theorem  3.4.5  fits 

within  the  union  of  these  two  excluded  regions.  Since  there  can  be  no  j 

separator  segment  through  p  p  ,  by  theorem  3.4.5,  edge  p  p  ,  the  shortest 

edge,  must  be  in  a  MWT. 
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Figure  21.  The  shortest  edge  is  in  a  MWT 


3.4  Delaunay  Triangulation  vs  MWT 


The  Greedy  Triangulation  has  been  shown  to  be  a  poor  approximation  to 

the  Minimum  Weight  Triangulation,  in  that  the  ratio  of  the  total  length  of 

the  Greedy  Triangulation  of  a  set  of  N  points  to  that  of  the  MWT  can  be  as 
1/3 

large  as  0(N  )  [5],  The  Delaunay  Triangulation  (DT)  is  the  straight-line 

dual  of  the  Voronoi  diagram  [1];  we  present  now  an  example  for  which  the 
ratio  of  the  length  of  the  DT  to  that  of  the  MWT  of  a  set  of  N  points  is 
roughly  N/2,  thus  showing  that  the  Delaunay  Triangulation  is  not  "approxi¬ 
mately  optimal". 

In  Figure  22,  five  points  are  closely  packed  on  an  arc  of  a  circle 
with  radius  r,  and  a  sixth  point  is  at  the  center  of  the  circle.  The 
Voronoi  diagram  and  Delanunay  Triangulation  are  shown  as  dashed  and  solid 
lines,  respectively.  If  the  five  points  were  so  closely  packed  that  the 
distance  between  the  furthest  two  were  negligible  as  compared  to  r,  then 
the  total  length  of  the  triangulation  would  be  roughly  5r.  In  general, 
for  N  points,  the  total  length  could  be  arbitrarily  close  to  (N-l)r. 


Figure  22.  Voronoi  diagram  (dashed  lines),  and 

Delaunay  triangulation  (solid  lines). 
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Figure  23  illustrates  the  Minimum  Weight  Triangulation  for  the  same 
set  of  points.  With  the  points  on  the  arc  closely  packed,  the  total 
length  can  be  arbitrarily  close  to  2r,  regardless  of  the  total  number  of 
points. 


Figure  23.  Minimum  Weight  Triangulation. 


Thus,  the  preceding  illustrates  a  construction  for  which  the  total 
length  ratio  between  DT  and  MWT  can  be  as  large  as  roughly  (N-l)/2,  for 
the  N^3,  and  the  Delaunay  Triangulation  need  not  be  a  good  approximation 
to  the  Minimum  Weight  Triangulation. 
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